EPJ manuscript No. 

(will be inserted by the editor) 



Anisotropic velocity distributions in 3D dissipative optical lattices 



J.Jersblad^, H.EUmann^, L.Sanchez-Palencia^, and A.Kastberg^''^ 
^f^' ^ Department of Physics, Stockholm University, S-106 91 Stockholm, Sweden 

o ' ' ' ■ - ' ■ ' - - ' - 

o 



0^ 



> 

o 

o 
m 
o 



X 



Laboratoire Kastler-Brossel, Departement de Physique de I'Ecole Normale Superieure, 24 rue Lhomond, F- 75231 Paris cedex 
05, France 

^ Department of Physics, Umea University, S-901 87, Umea, Sweden 

Received: date / Revised version: date 

Abstract. We present a direct measurement of velocity distributions in two dimensions by using an ab- 
sorption imaging technique in a 3D near resonant optical lattice. The results show a clear difference in the 
velocity distributions for the different directions. The experimental results are compared with a numerical 
3D semi-classical Monte-Carlo simulation. The numerical simulations are in good qualitative agreement 
with the experimental results. 

PACS. 32.80.Pj Optical cooling of atoms; trapping 



1 Introduction 



An optical lattice is a periodic optical light shift potential 
created by the interference of laser beams in which atoms 
can be trapped. Usually one distinguish between two types 
of lattices, near-resonance optical lattices (NROL) J] and 
far-off resonance lattices (FOROL). In the later type, an 
atom can only be trapped, whereas the former (the one 
considered in this paper) also exhibits an inherent cool- 
ing mechanism (Sisyphus cooling). The Sisyphus cooling 
mechanism in an NROL has been the subject of extensive 
research due to its high cooling efficiency, but also since 
an optical lattice is a very pure quantum system suitable 
for fundamental studies of atom-light interaction. 

Theoretical studies of the atomic motion in NROLs 
have been done in ID and 2D, both analytically and nu- 
merically. The extension to 3D configurations is however 
cumbersome. Analytical solutions become unwieldy and 
numerical simulations require long computation time, es- 
pecially for high angular momentum transitions. Thus very 
few detailed studies have been made in 3D. An exception 
is the work by Castin and M0lmer W who studied spatial 
and momentum localization via full quantum Monte Carlo 
wavefunction simulations in the case of optical molasses. 

Measurements of temperature have been made on 3D 
NROLs by our group j3 4 , and by groups at NIST and 
in Paris [H]. In all these experiments, and in this work, 
the kinetic temperature is derived from measured veloc- 
ity distributions along one axis and is defined as a direct 
measure of the kinetic energy through 



where 



is the mean square velocity of the released 



T 



atoms, fee is the Boltzmann constant and M is the atomic 
mass. A robust result in all studies is that the temperature 
scales linearly with the irradiance divided by the detun- 
ing, that is linearly with the light shift at the bottom of 
the optical potential {U). This is in excellent qualitative 
agreement with ID-theoretical predictions [J]. 

Nevertheless, in several works (for example [5] and 
this work) a four laser beam configuration results in a face 
centered tetragonal lattice that cannot simply be reduced 
to three ID cases. Indeed, all spatial directions are not 
equivalent (see section ITT|l and the particular geometry of 
the lattice has to be taken into account. In a recent paper 
by the Grynberg group .8 , the dependence of tempera- 
ture and spatial diffusion on geometric parameters con- 
trolling the lattice spatial periods (lattice constants) in 
different directions was studied. For different laser beam 
configurations producing the NROL, the temperature and 
spatial diffusion coefficient were measured for tetragonal 
lattices (see section|OJ with different aspect ratios, i.e. as 
a function of lattice constants. It was shown that the spa- 
tial diffusion coefficient strongly depends on the direction. 
The temperature, which was measured in one direction, 
was found to be independent of the lattice spacing. The 
difference between spatial directions lies not only in the 
lattice constants, but also in the modulations of the laser- 
atom interaction parameters (optical potentials and opti- 
cal pumping) in such a way that different behaviors of the 
temperature along different axes is possible. In [H| the Sisy- 
phus cooling effect in a 3D tetragonal NROL was studied 
theoretically. With a simplified choice of atomic angular 
momentum, it was shown by a semi-classical Monte-Carlo 
calculation that the temperature along a given coordinate 
axis is independent of the lattice constant, but indeed dif- 
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ferent along different directions. For the same geometry as 
considered here, the hnear scaling parameter of the tem- 
perature differs by a factor of 1.4. Moreover, a compari- 
son between |H] and H] suggests such an anisotropy of the 
velocity distribution. In both experiments, the direction 
of measurement coincided with the direction of gravity, 
but this direction did not correspond to the same lattice 
axis. It turns out that these works yield a quantitative dis- 
crepancy. The derived temperature was found to be linear 
with U with proportionality constants of 12 uK/E'r and 
24 uK/E'r (in and jH] respectively), where i?R is the 
recoil energy ^. The difference in scaling factor called out 
for a more thorough investigation, which would rule out 
any systematic error. 

This work aims at a direct comparison between the 
kinetic temperatures along different directions in a 3D 
NROL. Measurements of velocity distributions along dif- 
ferent directions were made for different lattice parameters 
(potential depth and detuning) by absorption imaging of 
an expanding atomic cloud. The experimental results are 
compared with a 3D semi-classical Monte-Carlo simula- 
tion performed for the actual atomic angular momentum. 

The paper is organized as follows. In section lTTI we de- 
scribe the experimental set-up. The experimental data is 
presented with derived kinetic temperatures in section 
In section |31 we describe the numerical calculation and 
present the result for the kinetic temperatures. In sec- 
tion 0] we discuss the results from the experiment and the 
simulations. Finally, in section [3 we draw conclusions on 
our work. 



2 Experiment 

2.1 Experimental setup 

Initially, a magneto-optical trap (MOT) is loaded with 
N 2 ■ 10® cesium atoms (^■^•^Cs) from a chirped deceler- 
ated atomic beam in 4 s. This gives a peak number den- 
sity of no w 5 • 10^° cm"3. The MOT operates at the 
(Fg = 4 Fe = 5) transition at 852 nm (the D2 hue), 
where F is the total angular momentum quantum num- 
ber. Due to off-resonant excitation to i^e = 4, a repumper 
beam resonant with the {Fg = 3 Fe = 4) transition is 
also used. After turning off the loading, the atoms are fur- 
ther cooled in an optical molasses for about 20 ms. From 
the optical molasses, an atomic cloud at a temperature of 
T = 3 /iK is loaded into the optical lattice with a transfer 
efHciency of about 50%. The filling factor of the lattice is 
around 0.2 %. The optical lattice beams are red detuned 
from the (Fg = 4 ^ Fe = 5) resonance, typically between 
Z\5 = -lOF and = -40F, where F/27r = 5.2 MHz is 
the natural linewidth. The atoms equilibrate in the lat- 
tice for 10 ms and are then released by turning off the 
optical lattice beams, with an acousto-optical modulator 
(AOM) in less than 1/is, followed by a measurement of 

^ The recoil energy E-^ = (ftfc)^/2M, where k = 2tv/\ is the 
wave vector, A is the wavelength of the light, and M is the 
atomic mass. 
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Fig. 1. Beam configuration of the 3D lin _L lin optical lattice. 
Two beam pairs in the xz- and yz-planes respectively, orthog- 
onally polarized along the y- and x-axes respectively, make an 
angle — 45° with the z-axis. 



the kinetic temperature. This short falltime of the AOM 
avoids adiabatic release of the atoms in the optical lattice. 

The optical lattice is a 3D generalization of the ID 
lin_Llin configuration created by two orthogonally polar- 
ized pairs of laser beams that propagate in the yz- and 
a;z-planes respectively . The angle between the beams of 
each pair is 90°, and each beam forms an angle of 6* = 45° 
with the (vertical) quantization (z-) axis (see figure^. 
This results in a tetragonal structure with alternating 
sites of pure cr"*"- and cr~-light, where potential minima 
are formed. From figure^ h is clear that directions x and 
y are equivalent but that direction z is different. It fol- 
lows that the optical pumping rates and the light shift 
modulations are different along z compared to x or y. In 
figure[21we plot the projection along x and z of the lowest 
adiabatic potential, which is where the atoms spend most 
of their time pp. Two main anisotropic properties arise. 
First, the lattice constants Uz = A/(2-\/2) and ax,y = X/V^ 
are different. Second, the shapes of the potentials are also 
clearly different. In particular, they show different poten- 
tial barriers to escape adiabatically from a potential well 
(lower along the z-direction than what it is along the x- 
and y-directions by a factor of 1.65) and show different 
reduced oscillating frequencies (wia^/A) at the bottom of 
the potential wells. 

The velocity distributions along the z- and a;-axes are 
measured using a well known absorption imaging tech- 
nique 10 . After release from the lattice, a short (50 /zs) 
resonant probe pulse (Fg = 4 ^ Fe = 5) hits the atomic 
cloud. The irradiance of the probe pulse is / <C /oi where 
Iq ~ 1.1 mW/cm^ is the saturation irradiance. The shadow 
in the probe beam is imaged onto a CCD camera. By cap- 
turing images at different time delays after turning off the 
optical lattice beams, we extract the different spatial den- 
sity distributions from which velocity distributions can be 
derived. The velocity distribution in the z-direction (di- 
rection of gravity) was compared to the results obtained 
with a " time-of-flight" (TOF) method JT] , showing good 
agreement. 
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Fig. 3. Typical 2D density profiles acquired at two different times r after releasing the atoms from the lattice. The left image 
shows an atomic cloud after r = 12.8 ms expansion together with density profiles in the z- and x-directions. The right image 
shows an atomic cloud after t = 36.8 ms. 




Fig. 2. Lowest adiabatic optical lightshift potential projected 
in the a;2:-plane in units of the optical wavelength, A. 



2.2 Measured Kinetic Temperatures 

The 2D projection (in the xz-plane) of the expanding 
cloud is recorded at two different time delays, ri_2, af- 
ter extinction of the optical lattice beams. Typical values 
are t\ = 12 ms and — 35 ms. Examples of 2D den- 
sity profiles are shown in figure |31 together with Gaus- 
sian fits to the spatial density profile along x and z. Ex- 
cellent agreement with Gaussian distributions is found. 
From the fits, we extract the rms radius, di, (i = a;, 2), 
of the clouds which increases with time, according to 

The kinetic temperature in different directions is de- 
fined as 



In figure 21 we plot derived kinetic temperatures along x 
and z for three different detunings, as a function of [/q, 
which is the modulation depth of the diabatic optical po- 
tential. Here, C/q is defined as 













1 + 


(1) 




2 




2ziiJ 



(3) 



where i?^ = {F'^ /2) / {I / If^) is the square of Rabi fre- 
quency and the irradiance is / = 8/boam (-^boam is the 
irradiance of a single beam), at the center of a potential 




Fig. 4. Tx (filled) and Tz (open) as a function of modula- 
tion depth, Uo/Er, for three different detunings [A^ = —IQF 
(squares) , —2Qr (triangles), —ZQF (circles)). The solid and 
dashed line are linear fits to the data. 



well. For sufficiently high irradiances and temperatures, 
it is obvious that the universal scaling with Uo prevails 
for each direction. However, this scaling with C/q is clearly 
different for different directions. For large ?7o, the temper- 
ature along z is found to be significantly smaller than the 
temperature along x. Linear fits to the data yield 

= (0.55 + 0.022(C/o/%)) (4) 
^ (0.62 + 0.012(C/o/%)) ^JK. (5) 

That is, the ratio between the scaling parameters along x 
and z is determined to be 1.8 (0.3). However, at low mod- 
ulation depths and low temperatures, the temperatures 
are found to be approximately the same along z and x. 

3 Numerical Simulations 

3.1 Theoretical framework 

We have performed semi-classical Monte-Carlo simulations 
in 3D for the actual (Fg = 4 — > i^e = 5) transition of ^'^^Cs. 
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The main features of the method have been discussed else- 
where so here we just recall the main elements and 
peculiarities for our multidimensional configuration. 

The optical Bloch equations (OBE), which describe 
the evolution of a sample of two-level atoms (with Zee- 
man degeneracy) coupled to both laser fields and vacuum 
modes, are the starting point of the analysis. Because of 
the cooling effects and the decoherence due to photon scat- 
tering, the atomic cloud dynamics can be reduced to a 
semi-classical picture for a large range of lattice parame- 
ters The OBE are therefore converted into a set of 
coupled semi-classical Fokker-Planck equations (FPE) via 
Wigner transforms. Projecting the FPE onto the position- 
dependent adiabatic states base |<?m(r)) (see appendix 
and neglecting the coherence terms which are unimpor- 
tant in a semi-classical description, one gets a new set of 
FPE only involving the local populations of the adiabatic 
states^. 

By physical interpretation of the FPE, it follows that 
the atomic cloud dynamics can be reduced to internal 
state transitions via optical pumping at a rate 7n,m from 
|<?„) to |<?m), and the evolution of each atom in a given in- 
ternal |<?m)-state due to deterministic forces. These forces 
are first of all due to the optical potential modulation (- 
VL^to) and secondly, due to the radiation pressure force 
(F). Moreover, the atomic cloud undergoes momentum 
diffusion due to photon scattering. 

It is then straightforward to show that the FPE solu- 
tion is formally equivalent to the integration of a set of 
Langevin equations interrupted by internal states quan- 
tum jumps, each one accounting for the random trajectory 
of a single atom. The quantum jumps are taken into ac- 
count by generating a random number r at each time step 
which is compared to the transition probability ^m,n 

dt 

from \<l>rn) to (with n ^ m) during the time step dt. 
In the following, we define r„ „ as 1 if a quantum jump oc- 
curs from n to m and otherwise. Between two quantum 
jumps, the elementary evolution of the atom is 



dRit)^^dt 



(6) 



dP (t) = -VU,ndt + r„,„ ((5p„ „ + r„,„,dt) 



)dt, (7) 

where R and P are the atomic position and momentum 
respectively. The Hamiltonian force, {-'VUm), is derived 
from the adiabatic potential in state |^m), and F„.m is 
the average radiation pressure in case of a quantum jump 
from n to m (if m = n, no jump occurs). The momentum 
diffusion is determined by random values: the momentum 
kick undergone by the atom in case of a quantum jump 
from n to m, <5p„ ^ and the recoil mean force in the ab- 
sence of a quantum jump, £„. Note that (5p„ „ and 



Note that the adiabatic approximation is justified by the 
fact that the adiabatic state splittings are generally greater 
than the motional couplings in the regime of deep potentials. 



are related to the position-dependent coefficients appear- 
ing within the FPE. In a (/z-indexed) space base where 
the momentum diffusion matrix {Dn^m} (see appendix IB|| 
is diagonal, the first two moments of tm and (5p„ „ read 



dt 



(8) 



(/„^,> - and iifiif) 
(^K,m> = and (((5p^:_,„ 

3.2 Numerical results 



The numerical simulations are performed for a typical 
sample of 300 independent atoms. For the lattice param- 
eters considered in this work, the kinetic energy reaches 
steady-state in a time of approximately 4000 /F , where 
r — rsQ/2 is the total scattering rate and sq is the satu- 
ration parameter (see appendix^l. The averages of the ki- 
netic energies in steady state in the x-, y- and z-directions 
provide the kinetic temperatures in the corresponding di- 
rections, 

T. ^ (9) 

The simulations were made for three different detun- 
ings (Z\5 = -lOr, -20r,-30r). For each detuning we 
acquired velocity distributions, in each direction, at six 
different modulation depths. Note that the chosen mod- 
ulation depths are much higher than in the experiment 
since the semi-classical model breaks down when the mo- 
mentum distribution becomes too narrow. This is because 
deep modulation depths are required to avoid non-adia- 
batic motional couplings between adiabatic sublevels that 
are not included in our treatment ^|. Moreover, the time 
to reach steady state increases for low modulation depths. 
However, the linear scaling should still hold. The results of 
the numerical simulations are shown in figure^! Here, the 
kinetic temperature is plotted as a function of modulation 
depth for the detunings mentioned above. The tempera- 
ture scales linearly with the light shift independently of 
the detuning according to 



T, cx OMbiUo/ER) nK 
Ty cx 0m5{Uo/EB,) /xK 
T, cx 0M3{Uo/Eb.) 



(10) 

(11) 
(12) 



As in the experiments, the results of the simulations show 
a clear difference in scaling of the kinetic temperature 
along the z-axis compared to the x- and y-axes, here, by 
a factor of 2.7. 



4 Discussion 

The results from the experimental work and the numerical 
simulations are compiled in tabled A comparison shows a 
quantitative excellent agreement between our experiments 
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Fig. 5. Kinetic temperature, along the x- (circles), y- 
(squares) and ^-direction (filled triangles) as a function of mod- 
ulation depth, Uq/Eyi, for three different detunings (Zis = 
— lO-T, —20r, —30r). The solid and dashed lines are linear fits 
to the data. 



Table 1. The scaling parameter ^x,y,z, (in units of nK/E-^), 
in the equation Tx,y,z ~ To + ^x,y,zUo for different studies. The 
experimental errors of the slope for this work is the quadratic 
sum of the statistical error and an estimated maximum sys- 
tematical error. The errors in the simulation is the statistical 
error from the fit. 

ref [3 ref this work this work 

(experimental) (simulations) 

^x - 24(2.4) 22(3.5) 35(1.2) 

- - - 35(1.2) 

12(1.2) - 12(2.5) 13(1.0) 



tice wells in agreement with former experimental inves- 
tigations, for instance "5^, and full quantum Monte-Carlo 
simulations 2 . The difference in the scaling factors is pro- 
portional to the difference in the modulation depth of the 
lowest adiabatic optical potential in the corresponding di- 
rections. Therefore we conclude that it is this difference 
which induce anisotropic kinetic temperatures in the op- 
tical lattice. This conclusion is not incompatible with the 
results of |S| in which the steady-state kinetic temperature 
was measured for different lattice constants showing that 
the steady-state kinetic temperature was independent of 
the lattice spacing, because the geometrical anisotropy in 
the lattice do not reduce to a simple scaling factor between 
directions x,y and z. 

At low modulation depths, the lattice reaches a min- 
imum temperature followed by a sharp increase in tem- 
perature, usually called decrochage. When laser cooling is 
still effective there exists a region where the temperature 
is isotropic. However, this region is difficult to analyze for 
several reasons. For instance, at low modulation depths 
the atomic localization in a trapping site is less strong, 
and thus the anharmonicity of the potential well becomes 
more important. This could lead to an increased coupling 
between the different motional directions and also a broad- 
ening of the vibrational levels, i.e. increasing the tunneling 
rate in the lattice. Another effect that must be taken into 
account at low modulation depths is increased spatial dif- 
fusion |8I16| . This means that the loss rate of the atoms 
in the lattice becomes larger, and thus the signal-to-noise 
in the absorption images decreases. Furthermore, if the 
thermal expansion of the atomic cloud in the recorded ab- 
sorption images is small compared to the size of the cloud, 
due to spatial diffusion, there will be large uncertainties 
in the extracted temperatures. 



and former studies in which the kinetic temperature was 
measured along a; or along z ■ The numerical simula- 
tions also reproduce the difference in scaling parameter for 
different directions was measured in the experiments, and 
confirms the appearance of a discrepancy between kinetic 
temperatures along x-y and z. 

The inherent cooling process in an optical lattice for 
atoms with kinetic energy Ek > Uq is Sisyphus cool- 
ing. This process was explained by Dalibard and Cohen- 
Tannoudji in |7] in the case of a theoretical transition 
{Jg = 1/2 ^ Je = 3/2). The Sisyphus cooling cycle occurs 
until the atomic kinetic energy is lower than the potential 
barrier in a particular direction and thus does not depend 
on any other anisotropy (the lattice spacings for example) . 

However, for higher angular momentum transitions, 
the cooling process does not stop because other relaxation 
processes than standard Sisyphus cooling could still occur 
|13I15| . For example, atoms in bound states within a lat- 
tice well can be excited to unbound states, followed by 
decay to lower lying vibrational states. We find that the 
atomic kinetic energy is £^k ~ and thus that the 

atoms are very well localized at the bottom of the lat- 



5 Conclusions 

We have measured the velocity distributions in a 3D op- 
tical lattice of cesium along two non-equivalent directions 
as a function of lightshift (t/o). In agreement with previ- 
ous works, the kinetic temperature scales linearly with Uq. 
As an original result, we have found that the distributions 
are clearly anisotropic (with Tx^y > Tz). The experimental 
results are in good agreement with a 3D numerical Monte- 
Carlo simulation and we conclude that it is the modulation 
depth of the adiabatic optical potential that determines 
the steady-state kinetic temperatures. The anisotropy in 
kinetic temperature is not paradoxical. In fact the "ki- 
netic temperature" here is defined as a simple measure 
of the atomic kinetic energy (see Eq. (^) and not as a 
thermodynamical temperature. This is because thermal- 
ization in Sisyphus cooling do not result from energy ex- 
change beetween particles via collisions, but from atom- 
photon interactions. Our result show that no thermody- 
namical temperature can be defined for Sisyphus cooled 
atomic samples because of the violation of the equipar- 
tition theorem |17| . Our results can give important clues 
for a full understanding of the cooling mechanism in an 
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optical lattice. Furthermore, knowledge about the veloc- 
ity distributions in all directions is important in precision 
experiments utilizing optical lattices. 
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A Optical Bloch equations and adiabatic 
states 

This appendix aims at introducing the adiabatic states for 
a general J J + 1 transition atomic sample. Consider 
an atom of dipole operator D ~ P(d+ + d~) , with 
being the raising and lowering components of D, and V 
the reduced dipole moment. This atom interacts with the 
laser field 



El (r,i)-:^e(r)e— 



(13) 



, where Eq is the amplitude of the electric field, wl is the 
laser frequency and e(r) is a vector describing the spatial 
varying profile of the laser polarization. The operators A 
and Bq represent the hermitian conjugates of the opti- 
cal pumping cycles (absorption of laser photons followed 
by emission of stimulated or spontaneous photons respec- 
tively), and are defined as 



A = 


d 


•e*(r) 




Bq^ 


d 







where 



e± 



■eir) 
■eg 

with (7 = 0, ± or q = x,y,z 



V2 



and eo = 



(14) 
(15) 

(16) 



are the circular basis vectors. After elimination of the 
excited state in the low saturation regime. 



So 



Z\2 + r2/4 



«1, 



(17) 



the atomic sample dynamics is governed by the OBE 
involving the projection of the density matrix onto the 
internal state including an Hamiltonian part: 



2M 2 ^ ' 



(18) 



plus a relaxation part. In the semi-classical limit, the 
position-dependent adiabatic states are defined as the eigen- 
states of the light-shift operator Ti^^A: 

h^A (r) (r)) = (r) |<?„, (r)) . (19) 

Note that in general (r)) and Um (r) cannot be cal- 
culated analytically. 



B Dynamics coefficients for the Langevin 
equation 

In this appendix, we give the general expressions for the 
dynamics coefficients involved in the FPE and Langevin 
equations for Sisyphus cooling in the low saturation and 
semi-classical regime. The transition rate from state |<?„) 
to state \^m) (for m ^ n) is 



q=±,a 



(20) 



The average radiation pressure term in the direction i 

(i= x,y,z) is 



Km - -f^n Im {^m\d,Bl\^n){<l^n\Bg\<P„ 

and the momentum diffusion matrix is 



(21) 



- E {{^.n\dlBlWn){<l>n\Bg\1>„ 



- {<P„^\d^BlWn){^n\djBq\$^) 
+ c.c) (22) 



q=±fi 



where Sa^p is the Kronecker symbol (1 when a = (3 and 
else) and i,j denotes the spatial directions {x,y^z). Note 
that for the sake of simplicity, the spontaneous emission 
pattern is simplified in a way that the photons are re- 
stricted to be emitted only along the x-, y- and z-axes. 
This approximation is justified because the kinetic energy 
is expected to be greater than the recoil energy |18| . 
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